Chiral Dynamics of Deeply Bound Pionic Atoms 
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We present and discuss a systematic calculation, based on two-loop chiral perturbation theory, of 
the pion-nuclear s-wave optical potential. A proper treatment of the explicit energy dependence of 
the off-shell pion self-energy together with (electromagnetic) gauge invariance of the Klein-Gordon 
equation turns out to be crucial. Accurate data for the binding energies and widths of the Is and 2p 
levels in pionic 205 Pb and 207 Pb are well reproduced, and the notorious "missing repulsion" in the 
pion-nuclear s-wave optical potential is accounted for. The connection with the in-medium change 
of the pion decay constant is clarified. 
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Recent accurate data on Is and 2p states of negatively 
charged pions bound to Pb nuclei |j[ have set new stan- 
dards and constraints for the detailed analysis of s-wave 
pion-nucleus interactions. This subject has a long history 
Pi [| d> Mj Hi @7 Hi? culminating in various attempts to 
understand the notorious "missing repulsion": the stan- 
dard ansatz for the s-wave pion-nucleus optical potential 
in terms of the empirical threshold 7riV-amplitudes times 
the proton and neutron densities p P , n , supplemented by 
sizeable double-scattering corrections, still misses the em- 
pirically required repulsive interaction by a large amount. 
On purely phenomenological grounds, this problem can 
be fixed by simply introducing a sufficiently large nega- 
tive real part, HcBq, in the p 2 term of the pion-nuclear 
optical potential || . The arbitrariness of this procedure 
is of course unsatisfactory, also in view of the fact that 
Re_Bo cannot be calculated microscopically with sufficient 
accuracy, and that the high precision measurements of 
the real part of the pion deuteron scattering length |j] 
suggests a very small net ReB by analogy. 

In the meantime, this issue has been revived in a vari- 
ety of ways. An interpretation of the "missing repulsion" 
in terms of a possible in- medium change of the pion decay 
constant, as suggested in ref. [jlOj, appears to be remark- 
ably successful 0, [L2| but has been debated Jl|, [□]• 
Clearly, this concept needs further justification. In the 
present paper we show that a key to the understanding of 
low-energy pion-nucleus s-wave interactions lies in its dis- 
tinct energy dependence imposed by chiral symmetry, in 
combination with the "accidental" approximate vanish- 
ing of the isospin-even threshold 7rA^-amplitude. Another 
important feature, generally ignored in previous analy- 
ses, is the systematic incorporation of gauge invariance 
at all places where the pion energy oj appears explicitly, 
when solving the Klein-Gordon equation in the presence 
of electromagnetic interactions. 

Our framework will be in-medium chiral perturbation 
theory at two- loop order . But before going into tech- 
nical details, the following simplified treatment may be 
useful to illustrate driving mechanisms. 



Consider a zero-momentum n~ interacting with nu- 
clear matter at low proton and neutron densities, p p and 
p n . The in- medium pion polarization operator to lead- 
ing order in the nucleon densities is expressed in terms 
of the isospin-even and isospin-odd off-shell ttN ampli- 
tudes T ± (w) as IL(u) = -T~(u) Sp - T+(u) p, with the 
isoscalar and isovector nucleon densities, p = p p + p n and 
Sp = pp — p n - The spontaneous and explicit breaking of 
chiral symmetry implies the following leading terms of 
those amplitudes fll5[]: 
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where f n = 92.4 MeV is the pion decay constant, 
TOtt = 139.57 MeV is the (charged) pion mass and ajy — 
45 ± 8 MeV jl7| is the pion- nucleon sigma-term. The em- 
pirical observation that T+(m 7r ) = (-0.04±0.09) fm ~ 
fist sets the constraint (3 ~ crjv/m 2 . 

Next, consider the pionic in-medium dispersion equa- 
tion at zero momentum, ui 2 — m 2 — TI(ll>) — . Introduce 
an (energy independent) equivalent optical potential U 
by to 2 = m 2 + 2 U with U <C m^, and expand around 
lo = TOtt- This gives 
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where the last approximate step introduces the wave 

an n- 1 
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ference between second and last step in (^ 
leading order. Inserting ([!]) and assuming dp <C p one 
finds 
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which is the expression proposed in ref. JlO| . It involves 
the driving Weinberg- Tomozawa term in T~(tu), but with 
the pion decay constant renormalized (f n — > f£{p)) to 
leading order in the density p, in accordance with the 
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corresponding in-medium change of the chiral quark con- 
densate < q q > . 

To the extent that U represents part of the (energy 
independent) s-wave optical potential commonly used in 
the phenomenological analysis of pionic atoms, at least 
part of the " missing repulsion" is thus given a physical in- 
terpretation in terms of the reduced in-medium /* in the 
denominator of (||). Of course, rather than constructing 
the potential U and following the steps leading to (|[||): 
one can directly solve the Klein-Gordon (KG) equation 
with the full energy dependence of the polarization oper- 
ator H(u>). This is the procedure systematically applied 
in this paper, with proper recognition of gauge invariance 
in the presence of the electromagnetic field. 

The KG equation with Coulomb potential V c (f) < 
and total pion self-energy, Ht ot (io, r) reads 

[(w - V c f + V 2 - m 2 - II tot (^ - V c , f )] $(f ) = . (4) 

The total polarization operator expressed in terms 
of local proton and neutron densities, n tot (a;,r) = 
Iltot ; Ppffl, Pn(f))-> can be split into its s-wave and p- 
wave parts 

Iltot (w; P P , Pn) = n(w) + AIIs(w; p p , p n ) + U p (lo; p p , p n ) , 

where we separate explicitly the phenomenological s- 
wave absorption term quadratic in densities, 
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parameterized as in ref. B. Here M stands for the nu- 
cleoli mass. We use ImBg = 0.063 m~ A and Re£>o = 
as our standard set and discuss variations of Rei?o later. 
For the p-wave part Hp (a;; p p , p n ) we use the traditional 
Kisslinger form with inclusion of short-range correlations 
and parameters as specified in ref. B (set A). 

Given the smallness of the isospin-even irN scatter- 
ing amplitude T + (lo), double scattering (Pauli-blocking) 
corrections in H(lo) are well known to be important [|2|. 
When those are included the " phenomenological" s-wave 
pion polarization operator becomes [Q : 

Ilphcn^/Op,/^) = -T~(tJ)5p-T+ s (u,p) p, (6) 

with 

T+(co,p) = T+H - ([T+H] 2 + 2 [T-H] 2 ) . (7) 

The local Fermi momentum kp{r) = [Sir 2 p(r)/2] 1 / 3 is 
rewritten in terms of the local density p(r). Taking the 
polarization operator (|^) at the on-shell pion energy lo = 



U(lo) = n pho n(w = m„; p p , p n ) , 



(8) 



and neutron density distributions p p (r) and p n {r) are 
given as two-parameter Fermi functions pj (r) = poj [1 + 
exp((r — Rj)/a,j)]~ . The central density po,j is nor- 
malized to the total number of protons and neutrons in 
the nucleus. The proton radii, R p , are extracted from 
the nuclear charge radii following from the analyses of 
muonic atoms [|l9|, taking into account the finite pro- 
ton size < r 2 >= 0.73 fm 2 : i? p [ 205 Pb] = 6.66 fm and 
i? p [ 207 Pb] = 6.67 fm. Since the charge radii have not 
been measured for the complete chain of Pb isotopes, 
we have interpolated linearly between two neighboring 
measured isotopes. The diffuseness coefficient is taken 
the same for 2 05,207p| 3 ^ a ^ _ q p or ^ e neu- 

tron radii we use values from the proton-neutron rms- 
radius difference as obtained in the Brueckner-Hartree- 
Fock calculations of ref.@: i?„[ 205 Pb] = 6.94 fm and 
i?„[ 207 Pb] = 6.97 fm. We assume a n = a p . The numeri- 
cal input is close to that in refs. (|[ 

Solutions of the wave equation ([|) for Pb isotopes 
with the energy independent (threshold) input (||) for 
the pion-nuclcar optical potential are shown in Fig. |lj by 
open circles. The filled circles in Fig. |l| are the results 
obtained with the polarization operator 



II(w) = n p hon(w;pp(r),p„(r)) , 



(9) 



in which we keep the explicit energy dependence as given 
by the driving terms (l|). The energy dependence ef- 
fects are evidently important, moving the calculated re- 
sults closer to the data. Indeed, with the gauge in- 
variant introduction of the electromagnetic interaction 
in the pion polarization operator (via the replacement 
lo —> lo — V c (r)) the off-shell pion-nucleon scattering am- 
plitudes are probed at energies lo — V c (r) > m^. This 
increases the repulsion in T~(lo) and disbalances the can- 
celation between the sigma-term cr^r and the range term 
-(3lo 2 in T+(lo), giving T+(lo - V c (r)) < 0. Omitting 
the replacement lo — > lo — V c (r) in II(a;), we would have 
lo < m-n-, and this would reduce the repulsion in T^(lo) 
and turn on attraction in T + (lo), thus leading in the 
wrong direction. Taking both the energy dependence and 
the proper gauge invariant substitution via n(cij — V c (r)) 
is therefore an essential ingredient. 

After these qualitative considerations we proceed now 
to the systematic calculation of the pion polarization op- 
erator using in-medium chiral perturbation theory. Here 
we extend the results of ref. [|l6] at the two-loop level by 
taking into account the explicit (off-shell) energy depen- 
dence. The polarization operator has the form: 

n(w) = no(w) + n ds H + n rcl (^) + n cor (w) . (10) 

The first term corresponds to the linear density approx- 
imation: 



we recover the traditional form of the (energy inde- 
pendent) s-wave optical potential ||. The proton 



n M 



~T-{lo)8 P -T+{lo)p. (11) 
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FIG. 1: Binding energies and widths of deeply bound pionic states in the isotopes 207 Pb (left figure) and 205 Pb (right figure). 
Diamonds show the experimental data from Uncertainties in the extraction of the Is level for 207 Pb are indicated by 
different choices of a control parameter R as specified in [jjj]. The results for the polarization operator (^|^|) are depicted by 
open and filled circles, respectively. Triangles show the results obtained with the chiral polarization operator (|fo|). 



The isospin-even off-shell vrA^-scattering amplitude at 
zero pion momentum can be written in the following form 
(for uj > m v ): 



crjv - (3oj 2 3g 2 A ml 3g 2 A Q 2 m^C 
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with 7 = (gATrfn/M) 2 + ln(2A/m„.) . The cut-off scale 
A = 737 MeV ~ 8/^ is chosen to reproduce the central 
empirical value of the on-shell scattering amplitude at 
threshold T~(m v ) = 1.85 ±0.09 fm g§. We neglect here 
small additional counter-term contributions proportional 
to the third order low-energy constants dj of ref . p3| . 



The next term in ( jl(i| ) corresponds to the important 
Ericson-Ericson double scattering correction ||J general- 
ized to isospin asymmetric nuclear matter and off-shell 
pions: 



where (3 = g\/^M - 2c 2 — 2c 3 , cat = -Acim 2 — 
9g\ml/MTTf 2 and T im = oj 2 Q/(8nf$) . The nucleon 
axial- vector coupling constant has the value gA = 1.27. 
We introduce the abbreviation Q = ^/u> 2 — m 2 . The where k 
second-order low-energy constants 01,2,3 (f° r notations 
see ref. are tuned to the empirical values of the 

sigma-term [[l7|, — 45 MeV, and the ttN scattering 
length, T+(m 7r ) = 0. 

The parameter £ reflects freedom in the choice of 
the interpolating pion field in the effective chiral La- 
grangian jll| |2^]. It enters all interaction vertices with 
three and more pions. The one-loop correction to the 
(off-shell) pion self-energy in vacuum depends also on 
this parameter £. By requiring that the residue at the 
pion pole remains equal to one [pi] as it is implicit in 
the form of the KG equation (^) , one gets the constraint 

C = o. 

The isospin-odd off-shell 7r./V-amplitude at zero mo- 
mentum reads: 
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Q) + 2L(k p ,k n ,Q)}, (13) 



,-,.„ — (3tt 2 pp^) 1 / 3 refer to the proton and neu- 
tron Fermi momenta and 



L(kp,k n ,Q) — 4:kpk n (Q + 3fc + 3fc n ) 
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The third term in ( |iTj| ) is a small relativistic correction 
from the particle-hole (Born) diagram evaluated at zero 
pion momentum: 



n rc i(w) = 
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M~k 5 n ). (15) 



The last term in (11) represents the effect induced by 7r7r- 
interactions with two virtual pions being absorbed on the 
nucleons in the Fermi-sea, and by an additional two-loop 
correction |l(| |2^] : 
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FIG. 2: Dependence of energies and widths of pionic levels for 
205 Pb on parameters of the optical potential. Solid lines cor- 
respond to the variation of the neutron radius in the interval 
SR n — —0.2-^0.2 fm. Dashed lines show the variation with the 
pion nucleon sigma term (Tjv(0) = 20 -f- 65 MeV. Dotted lines 
represent the variation of Re£>o = (—0.5 ~ 0.5) IrrrBo. Arrows 
indicate the directions in which results move as the varied pa- 
rameters increase. Triangles and experimental points are the 
same as in Fig. |lj 
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The function H(k p ,k n ) consists of the last four terms 
written in eq. (12) of ref. ||(|. In the actual calculation 
the contribution H coi (uj) turns out to be negligibly small. 

Solutions of the wave equation (||) for Is and 2p levels 
using ( |To| ) are shown in Fig. ^ by triangles. They agree 
well with experimental data. We also find that the energy 



dep< indent polarization operator (|6 , 7 ) gives equally good 



results as (|10|) if the amplitude T~(uj) in (|l|) is extended 
to include the to 3 term in (|l2|), with the parameter 7 
tuned to the empirical value of the scattering length. 

In Fig. ^ we examine the dependence of our results for 
205 Pb on the less constrained parameters of the model. 
Variations of the neutron radius R n affect mainly the 
binding energy, whereas the sigma term ctn and Re£?o 
have a stronger impact on the level width. Note the 
strong correlation in the effects induced by changes of 
RcSo and <tn ■ 

In summary, we have demonstrated that the long 
standing issue of the "missing repulsion" in the s-wave 
pion nucleus potential can be naturally resolved by tak- 
ing into account the explicit energy dependence of the 
pion self-energy and the gauge invariant incorporation 
of electromagnetic interactions. The experimental data 
for Is and 2p levels in Pb isotopes are well reproduced. 
We have also clarified that, to leading order, the energy 
dependence effects can be interpreted in terms of the in- 



medium reduction of the pion decay constant, or equiv- 
alently, the renormalization of the intrinsic pion wave 
function in matter. 
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